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Aims. The production of saturated organic molecules in hot cores and corinos is not well understood. The standard approach is 
to assume that, as temperatures heat up during star formation, methanol and other species evaporate from grain surfaces and 
undergo a warm gas-phase chemistry at 100 K or greater to produce species such as methyl formate, dimethyl ether, and others. 
But a series of laboratory results shows that protonated ions, typical precursors to final products in ion-molecule schemes, tend 
to fragment upon dissociative recombination with electrons rather than just ejecting a hydrogen atom. Moreover, the specific 
f****) , proposed reaction to produce protonated methyl formate is now known not to occur at all. 

Methods. We utilize a gas-grain chemical network to probe the chemistry of the relatively ignored stage of hot core evolution 
during which the protostar switches on and the temperature of the surrounding gas and dust rises from 10 K to over 100 
K. During this stage, surface chemistry involving heavy radicals becomes more important as surface hydrogen atoms tend to 
evaporate rather than react. 

Results. Our results show that complex species such as methyl formate, formic acid, and dimethyl ether can be produced in 
large abundance during the protostellar switch-on phase, but that both grain-surface and gas-phase processes help to produce 
most species. The longer the timescale for protostellar switch-on, the more important the surface processes. 
Conclusions, xxx 
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1. Introduction chemistry. Various hot core mode ls llMillar et all 1 1991: 

, , , , , , „ ICharnlev et"al]ll992i ICaselli et al.l fl993t ICharnlev et alJ 

Various large molecules including methyl formate TnocITT iT 4.1,4.1 u 11 

1995J) have shown that large oxygen- bearing molecules 



(HCOOCH3) 

hot moleculai 


have been detected in a number of 
cores and corinos (IBlake et alJ 1987: 


Hatchcll et al. 


1998; 


Nummelin et alJ|2000: Cazaux et alJ 


2003: Bottinclli et al. 2004). These arc chemically rich re- 



such as methyl formate may be formed during the "hot" 

stage when large amounts of methanol are released from 

the grain mantles. Methanol is easily protonated by 

ions, and the resulting protonated methanol was thought 
gions in dense interstellar clouds which are warmed by , , r , 1 , , , , , , , ,1 , 

. ...... c to react with formaldehyde to produce protonated methyl 

an associated protostar and have typical densities of 10 r , , , , . , 

a . . , formate. Subsequent dissociative recombination with elec- 

- 10 5 cm 3 and temperatures on the order ot 100 K. , , , 1+4.1, t 4.1, • 

r trons would then complete the process. In this way, 

Observations suggest methyl formate abundances can be, 4. r 4.1. 1 r 4. uu jj 

°° J large amounts ol methyl formate could be produced on 

as high as ~10 x n(H2j in these regions. ,. , r m4 -in5 a. 4.1. 4.- r 4.1. 

° . . timescales ol 10 - 10 yr after the evaporation of the 

The standard view of the complex chemistry that per- 
tains in hot cores requires that dust grains build up icy 

mantles at early times, when the temperatur e is low (~10 . 

T r\ j ,1 • • r ii I7T~ — T — Tl rather more contentious. Recent quantum chemical calcu- 

K) a nd the core is m a state ol collapse (Brown ctaJJ , , . , , , r— — — — rfrrTTT, , 

r-T*— 1 t j , n e . ... lations carried out by Horn et al. ( 2004) have shown that 

1988). Later, the formation ol a nearby protostar quickly „ „ .„ ,7 r, TT . 

. , • the hrst stage ol the process - reaction between H 2 CO and 

the gas and dust, re-injecting the gram man- nu _ u+ „«,„,_+ „ + j„^„„ *_„+„j 



ices. 

However, this gas-phase mechanism has now become 



warms 

tie material into the gas phase, and stimulating a rich 



CH 3 OH^ - is highly inefficient at producing protonated 
methyl formate. This is due to a large potential barrier 
Send offprint requests to: R. T. Garrod, e-mail: between isomeric complexes, resulting in an activation en- 
rgarrod@mps.ohio-state.edu ergy barrier for the reaction on the order of ~15,000 K. 
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Horn et al. suggested other gas-phase routes to produce 
protonated methyl formate as well as two other complexes, 
adopting moderate efficiencies for each to recombine with 
electrons to produce HCOOCH3. Their chemical models 
(at 100 K) showed that even with optimistic rates, these 
reactions were unable to reproduce observed abundances. 

This considerable difficulty is compounded by increas- 
ing experimental evidence regarding the efficiencies of 
the dissociative recombina tions of large saturated ions. 
iGeppert et al.1 (|2005l l2006(l have investigated the branch- 
ing ratios of the recombination of protonated methanol, 
finding that the fraction of recombinations producing 
methanol is considerably lower than the 5 - 100 % typi- 
cally a ssumed in chemical m odels (see e.g. iLe Teuff et al.1 
|2000|) . IGeppert et all l|2006|) find that just 3 ± 2 % of re- 
combinations result in methanol, with more than 80 % 
resulting in separate C and O groups. This work is in 
agreement with many previous storage-ring experiments 
on the products of the dissociative recombination reac- 
tions involving other H-rich species, which are dominated 
by three-body channels. The strong implication for such 
larger molecules as methyl formate and dimethyl ether, 
which were hitherto expected to form from their proto- 
nated ions, is that their production via recombination 
must be, at best, similarly inefficient. Therefore, the gas- 
phase routes to such molecules must be considered inca- 
pable of reproducing astronomically observed abundances 
at temperatures of ^100 K. 

Numerous allusi on s have been made (e.g. 
Charnlev et all Il995t lEhrenfreund fc Charnlevl !000; 



Table 1. Desorption energies of selected species 



Cazaux et al.ll2003l IPeeters et al.1 2006j) to the possibility 



that methyl formate is produced on dust-grain surfaces 
at some point in the evolution of a hot core when dust 
temperatures are relati vely low. The reaction set of 
lAllen k Robinson! l|l977j) . a list of exothermic and plausi- 
ble grain-surface reactions, includes one such association 
between the surface radicals HCO and CH3O, producing 
HCOOCH3. Many of the reactions listed in this set have 
been incorporated into the gas-g rain chemical model o f 
The Ohio State University (OSU) l|Hasegawa et al.lll992|) . 
but only very few of those involve more than one large 
radical. This is because at the low temperatures of dense 
interstellar clouds (^10 K), such large radicals would be 
effectively stationary on the grain surfaces, prohibiting 
their reaction with anything but small, relatively mobile 
species like atomic hydrogen. At the high temperatures 
encountered in hot molecular cores, these radicals would 
no longer be present on grain surfaces at all, having 
thermally evaporated. However, in the light of the failure 
of gas-phase chemistry to reproduce methyl formate 
observations, we must ask how these radicals should 
behave at intermediate temperatures. Clearly such a 
physical state must exist in the timeline of hot core 
evolution, but is it long enough to significantly influence 
the behaviour of the chemistry? Since most chemical 
models involve only a rudimentary treatment of grain 
surfaces, and implement only a gas-phase chemistry, few 
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CO 


1150 


CH4 


1300 


HCO 


1600 


H 2 CO 


2050 


CH 3 


5080 


CH 3 OH 


5530 


HCOOH 


5570 


HCOOCH3 


6300 


CH3OCH3 


3150 



investigations of the warm-up phase of hot cores have 
be en carried out. 

IViti fc Williams! lll999l) explored the effects of the se- 
lective, time-dependent re-injection of grain mantle-borne 
species, using the same thermal evapor ation mechanism 
as em ployed in the OSU gas-grain code ijHasegawa et alJ 
119921) . They used a linear growth in the temperature of the 
hot core, starting from 10 K, over periods corresponding 
to the switch-on times for hydroge n burning of variously 
sized protostars, as determined bv iBernasconi fc Maederl 
They investigated the effect of selective re- 
injection of mantle species on gas-phase che mical evolu- 
tion. The follow-up work of lViti et al.l l)2004j) used results 
from temperature programmed desorption (TPD) labora- 
tory experiments to determine temperatures at which par- 
ticular species should desorb, taking into account phase 
changes in the icy mantles. Whilst these investigations 
shed new light on the dependence of the chemistry on 
(currently ill-defined) physical conditions, they necessar- 
ily ignored the effects of the grain-surface chemistry prior 
to re-injection. 

In this study, we use the OSU gas-grain code to 
model the warm-up phase of hot-core evolution, with the 
primary aim of reproducing observed HCOOCH3 abun- 
da nces. The code was previously used to model hot cores 
bv lCaselli et alJ l|l993j) . requiring constant grain tempera- 
tures. The code is now capable of dealing with the variable 
temperatures we expect in the warm-up phase. The grain 
chemistry allows for time-dependent accretion onto and 
evaporation (thermal and cosmic ray-induced) from the 
grain surfaces. Surface-bound species may react together 
- the majority of reactions are hydrogenations but reac- 
tions may take place between other atoms, and smaller 
radicals such as OH and CO, at rates determined by the 
grain te mperature and diffusion barriers. Following the 
work of iRuffle fc HerbstJ ll2001al) we allow surface-bound 
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molecules to be photodissociated by the (heavily extin- 
guished) interst ellar radiation field, and by the cosmic-ray 
induced field of IPrasad fc Tarafdarl l|l983^ . 

To this b asic chemical mod e l we add three new re- 
actions from lAllen fc Robinson] l|l977|) . allowing for the 
formation of methyl formate, dimethyl ether (CH3OCH3) 
and formic acid (HCOOH) from grain-surface radicals. We 
employ a new set of desorption and diffusion barriers rep- 
resenting water ice. We also investigate the e ffects of the 
warm- up on the gas-phase reactions set out in lHorn et alJ 
(2004J). Those reactions show a strong inverse temperature 
dependence, implying much larger rates at temperatures 
less than 100 K. 



2. Model 

2.1. Physical Model 

We adopt a two-stage hot-core model: in stage 1, the gas 
collapses from a comparatively diffuse state of nn — 3 x 
10 3 cm -3 , according to the isothermal free -fall collapse 
mec hanism outlined in lRawlings et all l|2002|) and derived 
from lSnitzerl l(l978l) . This process takes just under 10 6 yr 
to reach the final density of 10 7 cm~ 3 . 

In stage 2 the collapse is halted, and a gradual increase 
in dust and gas temperatures begins. W e base this phase o f 
our physical model on the approach of lViti et alJ l|2004j) . 
T hey used the observed protostellar luminosity function 
of iMolinari et al.1 l|2000l) to derive effective temperatures 
throughout the accretion phase of a central protostar, up 
until it reaches the zero-age main sequence (ZAMS), at 
which hydrogen burning takes place. They approximated 
the effective temperature profile as a power law with re- 
spect to the age of the protostar, then assumed that the 
temperature of the nearby hot core material should follow 
the same functional form. The contra ction times deter- 
mined bv iBernasconi fc Maederl l|l99rJ) were used to pro- 
vide a timescale for this process, on the order of 10 4 to 
10 6 yr for stars of 60 to 5 solar masses, respectively. 

Here, we adopt power law temperature profiles for the 
hot core material, with a power index of either 1 (lin- 
ear time-dependence) or 2. We allow the temperature 
to increase from 10 - 200 K over three representative 
timescales: 5 x 10 4 , 2 x 10 5 and 1 x 10 6 yr, approximat- 
ing high-mass, intermediate, and low-mass star-formation, 
respectively. The high gas density of 10 7 cm~ 3 implies a 
strong collisional coupling between the dust and gas tem- 
peratures, hence we assign them identical values at all 
times. 



1992). Similarly, desorption energies determine the rate of 
thermal and cosmic ray h eating-induced desorption (see 
IHasegawa & Herbst 1993). Diffusion barriers are set as a 
fixed fraction of the desorption energy of each species, as in 
previ ous work fe.g. IHasegawa et alll992!|Ruffie fc Herbst! 
2000). Both the diffusion rates and thermal desorption 
(evaporation) rates are exponentially dependent on the 
grain temperature, according to a Boltzmann factor. 

Of the desorption energies that have been used in pre- 
vious applications of the gas-grain code, some (including 
those for C, N, O, and S) correspond to a water ice surface, 
but most are given for carbonaceous or silicate surfaces. 
Here we construct a set of values corresponding to wa- 
ter ice fo r all the surface species in the gas-grain code, 
following ICuppen fc HerbsT l|2006j) . These authors anal- 
ysed various experimental data to obtain best values for 
H, H 2 , O2 and H 2 on non-porous ice. We use these val- 
ues and add to them data provided by Dr. M. P. Collings 
(private communica tion), which correspo nd to the TPD 
results obtained by ICollings et alJ l|2004j) . This gives us 
values for CO, C0 2 , CH 3 OH, C 2 H 2 , HCOOH, N 2 , NH 3 , 
CH 3 CN, SO, S0 2 , H 2 S, and OCS. (We ignore the value 
for CH4, which was su bject to experimental problems - 
see lCollings et al1l2004|) . We additionally assume that the 
OH value is half that for H 2 0, to take account of hy- 
drogen bonding. We use these skeletal values (along with 
unchanged atomic values, except for H) and extrapolate 
to other species by addition or subtraction of energies 
from members of the set, making the assumption that 
the energies are the sum of the fragments that compose 
the molecules. The net result is that most species possess 
desorption barriers around 1.5 times larger than in our 
previous models, whilst those species that should exhibit 
hydrogen bonding to the surface can have barriers as much 
as five times as large. Tabled shows a selection of impor- 
tant surface species and their desorption energies, Ed- A 
complete table of energy parameters is available from the 
authors on request. 

In previous models, diffusion barriers were set to ei- 
ther 0.3 or 0.77 of the desorption energies for most species 
(not a tomic or molecular hydrogen); see lRuffle fc HerbstI 
(2000). Here we adopt a blanket value for all species of 
0.5 Ed, based on an assortment of data on icy surfaces 
and our best estimate of the role of surface roughness and 
porosity. 

The reaction se t we em ploy derives from that used 
bv lRuffle fc HerbstI I 2001al) . which includes grain-surface 
photodissociation. To that set we add the follow ing three 
surface reactions from I Allen fc Robinson! 1 1977jl : 



2.2. Surface Chemistry 

For the grain-surface chemistry, we employ a new set of 
diffusion energy barriers and desorption energies for each 
species. Because the diffusion barriers determine the rate 
at which a species can hop between adjacent binding sites 
on the surface, the sum of these rates for two reaction 
partners determines the reaction rate I Hasegawa et al.1 



HCO + CH3O -> HCOOCH3, 
CH 3 + CH3O -> CH3OCH3, 
HCO + OH HCOOH. 



(1) 
(2) 
(3) 



These reactions we deem to be the most likely sur- 
face routes to each of the three molecules. During the cold 
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Table 2. Reactions from lHorn et alJ l|2004|) . with adjusted rates a 



Reaction 




fc(300 K) 


T-dependence, a 


fc(10 K) 






(cm 3 s" 1 ) 




(cm 3 s" 1 ) 


(4) CH3OH+ + H2CO 


-> CH3OH2OCH+ + hv 


3.1 x 10" 15 


-3 


8.4 x 10 -11 


(5) H 2 COH+ + H 2 CO 


-> H2COHOCH+ + hv 


8.1 x 1(T 15 


-3 


2.2 x 10~ 10 


(6) CH+ + HCOOH - 


> HC(OH)OCH+ + fei/ 


1.0 x 10 -11 


-1.5 


1.6 x 10" 9 



a k{T) = fc(300 K) x (T/300) 



Table 3. Initial abundances of H 2 and elements. 



Species i 


Ui/UH f 


H 2 


0.5 


He 


0.09 


C 


1.4(-4) 


N 


7.5(-5) 





3.2(-4) 


s 


1.5(-6) 


Na 


2.0(-8) 


Mg 


2.55(-6) 


Si 


1.95(-6) 


P 


2.3(-8) 


CI 


1.4(-8) 


Fe 


7.4(-7) 



+ a(b) = axlO b 



phase, the radicals HCO and CH 3 are formed by suc- 
cessive hydrogenation of CO, which mostly comes directly 
from accretion from the gas phase. The reactions H + CO 
— > HCO and H + H2CO — > CH3O each assume an acti- 
vation energy of 2500 K (|Ruffle fe Herbst! 1200 la|) . Large 
reserves of H2CO and CH3OH may build up through the 
hydrogenation process. However, at warmer temperatures 
such processes should be less efficient, due to the faster 
thermal evaporation of atomic hydrogen. Under these con- 
ditions, processes which break down H2CO and CH3OH 
may become important. 



2.3. Modified Rates 

Applications of the gas-grain code to cold dark clouds 
in previous publications have made use of "modified" 
rates in certain situations; specifically, problems arose 
when solving for species with average grain-surface abun - 
dances on the order of unity or below l|Caselh et alJll998h . 
Compa rison with results from simple Monte Carlo tech- 
niques llTielens fc Hagenlll982ft indicated that atomic hy- 
drogen in particular reacted too quickly when using the 
rate-equation method employed in the gas-grain code. In 
effect, the hydrogen was reacting more quickly than it was 
actually landing on the grain surfaces. This problem was 
fixed by limiting the rate of diffusion of the hydrogen to no 



more than the rate of accretion. Additionally, when in the 
"evaporation limit" (i.e. when the rate of evaporation ex- 
ceeds the rate of accretion), the diffusion rate was limited 
to the rate of H evaporation instead, since in this case it is 
the evaporation rate which is the main determinant of sur- 
face hydrogen abundance. Whilst IShalabiea et alJ (^998) 
applied the modified rate technique to all species, the ef- 
fect was found to be small except for atomic H, at the 
lo w temperatures of dark clouds (~10 K). The analysis 
bv lKatz et all l)l999h of earlier experimental work subse- 
quently showed that surface hydrogen hops between bind- 
ing sites thermally, rather than via quantum tunnelling, 
indicating a much reduced diffus ion rate for atomic hy- 
drogen . The use of this result by iRuffie fe Herbstl (|2000L 
l2001al) showed that the need for modified rates even for 
atomic hydrogen was marginal. 

However, in this work, we deal with much higher tem- 
peratures - at which even the heaviest species can become 
highly mobile. This means that any species could react 
at a faster rate than it is accreted or evaporated when 
some critical temperature is reached. But this does not 
necessarily mean that modified rates should be invoked 
for all species and reactions. Many of the heavy radicals 
in which we are interested here are formed mainly by sur- 
face reactions, and are therefore not governed by accretion 
in the way that many species are whilst at low tempera- 
tures. Also, some species can be produced by cosmic ray 
photodissociation of an abundant surface molecule. These 
complications suggest that any need to apply the modified 
rates is highly specific, both to individual species and to 
particular brief periods for those species. Added to this 
is the fact that the modified rates are untested against 
Monte Carlo codes at temperatures much greater than 20 
K, meaning that there is no certain need to modify the 
rates at all in those cases. 

Hence, we modify only reactions involving atomic hy- 
drogen. It is the only species that requires modification 
at less than 20 K, and is easily evaporated above 20 K, 
making it less important in the surface chemistry at such 
temperatures. As a result, the modified rates mainly come 
into play in the low-temperature collapse phase. The rates 
of surface reactions involving at omic hydrogen are mod - 
ified in the manner outlined bv IRuffie fc Herbstl l)2000h . 
Additionally, we modify the rates of hydrogen reactions 
t hat have activa t ion en ergy barriers, following the method 
of lAikawa et alJ 1 2005|) . This addition is necessary because 
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at the high densities involved, the quantities of reactant 
that build up can become large enough that, even with the 
presence of the activation energy barrier, the rate at which 
hydrogen reacts is greater than the rate at which it is be- 
ing deposited on the surfaces. To avoid such a run-away 
hydrogen chemistry, if the product of the activation bar- 
rier term, the surface abundance of the (non-hydrogenic) 
reactant, and the diffusion rate, is greater than the faster 
of the rates of accretion and evaporation of atomic hydro- 
gen, then that product is replaced in the rate equation 
with the accretion o r evaporation rate. For a full discus- 
sion of this issue see lCaselli et al.l (|2002). 



2.4. Gas-Phase Chemistry 



Table El shows reactions previously included in lHorn et alJ 
( 2004) that lead, without activation energy, to protonated 
methyl formate and to two other methyl formate precur- 
sors. Two of the ionic products require significant struc- 
tural re-arrangement in a subsequent dissociative recombi- 
nation reaction to produce methyl formate. The rate coef- 
ficient given to the first reaction in the table (labelled 4) is 
unchanged; to the other two (5 and 6) we attribute some- 
what weaker temperature dependences due either to com- 
petition or a better estimate. At a temperature of 10 K, 
all three reactions are significantly faster than at the 100 
K at which they were tested by Horn et al.; reaction (6) 
reaches approximately the collisional (i.e. maximum) rate 
coefficient at this temperature. Hence, we expect contri- 
butions from these reactions to be strong during the early 
stages of the warm-up phase, if the abundances of the re- 
actants are significant. We allow the products of reaction 
(6) to recombine with electrons to form methyl formate in 
5 % of reactions, two times lower than assumed by Horn 
et al. The protonated products of reactions (4) and (5) 
are allowed to form methyl formate in just 1 % of recom- 



binations, because each of these channels would require a 
more drastic structural change within the complex. 



We also include appro ximately 200 gas-pha se reactions 
from the UMIST ratefile l|Le Teuff et al.l2000|) which were 
not present in the (predominantly low temperature) OSU 
reaction set. These reactions have activation energies, but 
the greater temperatures seen in hot cores justify their 
inclusion here. 



The initial gas-phase matter is assumed to consist of 
atoms with the exception of molecular hydrogen. Shown 
in t able |3] these values derive from those selected by 
IWakelam fc Herbst! |2006) from the most recent observa- 
tional elemental determinations for diffuse clouds. Since 
the high temperatures attained in hot cores are capable 
of driving off the entirety of the granular icy mantles, we 
choose relatively undepleted (i.e. diffuse cloud) values for 
the lighter elements, in contrast with prev ious applica- 
tions of the gas-grain code (most pertinently. ICaselli et alJ 
Il993l) . We allow heavier elements to be depleted by an 
order of magnitude, on the assumption that the remain- 
der is bound within the grain nuclei. This also ensures 
that frac tional ionizatio ns stay at appropriate levels, de- 
fined bv iMcKeel lll989h as X(e) ~ 10- 5 n{R 2 )~ 1/2 ■ The 
behaviour of sulphur on grain surfaces is not fully un- 
derstood, and the form it takes on/in the grains is not 
known; observations have so far failed to detect H2S on 
grains, and some authors suggest P CS may be the domi- 
nant form l|van der Tak et al Jl2003^ . Moreover, the levels 
of sulphur-bearing species detected in the gas phase of hot 
cores are typically n ot as high as the levels detected in dif- 
fuse clouds (see e.g. lWakelam et al"ll2004a|) . We therefore 
de plete sulphur by an order of magnitude from the values 
of lWakelam fc Herbst! fcOOd) . 
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Fig. 2. Temperature profiles for a protostellar switch-on 
time of t — 2 x 10 5 yr, adopting linear (Ti) and t 2 time 
dependences (T2). 



phase. Once formed, CO2 does not take part in further 
reactions, but may be photodissociated. 

3.2. Stage 2: Warm-up Phase 

We first explore the behaviour of a "standard" physical 
and chemical model, before discussing other variations and 
free parameters. Figure shows two temperature profiles 
used in the warm-up phase of the model; our "standard" 
model uses the AT2 profile with a protostellar switch-on 
time of 2 x 10 5 yr, corresponding to the final time acheived 
in the model, tf. All of reactions (1) - (6) are enabled. 
Figure [21 shows plots of pertinent gas-phase and grain- 
surface species over the final "decade" of evolution; i.e. 
2 x f 4 - 2 x f 5 yr. Over this period, temperatures begin 
to rise significantly above 10 K, and ultimately reach the 
final value of 200 K. 



3. Results 

3.1. Stage 1: Collapse 

All models investigated here use the same stage 1 collapse 
model, hence the initial conditions for the warm-up phase 
are the same in all cases. Figure ^ shows the evolution of 
some important gas-phase and grain-surface species, with 
respect to n#, during the collapse from 3 x 10 3 - 10 7 cm~ 3 . 

Initially, the collapse is slow, and most of the signifi- 
cant increase in density takes place during the final 5 x 10 5 
yr. Over this period, we see that gas-phase reactions be- 
come much faster due to the greater density, whilst accre- 
tion rates also increase, allowing for an accelerated sur- 
face chemistry. Large amounts of H 2 0, CO, CH 4 , H 2 CO, 
CH3OH, H2S and NH3 have built up on the grain surfaces 
by the end of the collapse. Besides water ice, NH3 and CH4 
make up a large and fairly constant proportion of the ice. 
CH4 ranges from ~10 - 24 % of total ice throughout the 
collapse phase; NH3 is steady at around 15 %. Water, 
methane and ammonia ices are formed by successive hy- 
drogenation of O, C and N respectively. Carbon monoxide 
is predominantly formed in the gas phase and then ac- 
creted onto surfaces. At the very end of the collapse, CO 
deposition outpaces H 2 formation, due to the shortened 
supply of atomic oxygen. The model suggests that the 
outer 50 or so monolayers of ice should have a COiB^O 
ratio of approximately 6:4. Formaldehyde, which is formed 
by hydrogenation of CO, should comprise around 5-10 
% of these outer layers. 

The calculated abundances of ices at the end of the 
collapse are in reasonable agre ement with IR observa- 
tions fe.g. lNummelin et alJl200l|) except for CO2, which is 
greatly underproduced at 10 K. Solutions to th is problem 
have been advanced bv lRuffle fc HerbsH l|2001bh . Observed 
CO2 percentages (in dark cloud ices) are typically the 
same as or a little less than CO levels. Hence, even with 
efficient CO2 production, there should not be a drastic ef- 
fect on CO, H 2 CO, or CH 3 OH ice abundances in the cold 



3.2.1. Methyl Formate, Formic Acid, and Dimethyl 
Ether 

Methyl formate is formed in significant quantities, both 
in the gas phase and on grain surfaces. Figure shows 
the gas-phase abundance at late times - a double- 
peaked structure is apparent. The first, lesser, peak comes 
about solely from gas-phase formation via reaction (5). 
Formaldehyde injected into the gas phase from the grain 
surfaces at this time is protonated by reaction with H3 + 
and HCO + ions. At t — 8 x 10 4 yr, a temperature of 
~40 K obtains, giving a rate coefficient for reaction (5) 
of 3.4 x 10 -12 cm 3 s _1 . This is more than an order of 
magn itude larger than the rate at 100 K (cf. iHorn et alJ 
2004), a temperature typical of hot core chemical models. 
These factors combine to produce just under 10 -8 with re- 
spect to H2 of methyl formate; however, this abundance of 
HCOOCH3 does not manifest itself fully in the gas phase, 
since most of it is quickly accreted onto the grains. (The 
high density produces an accretion timescale of ^2000 
yr). This causes the sharp fall-off in gas-phase abundance. 
Reactions (4) and (6) (see table|5J) do not contribute signif- 
icantly to methyl formate production at any stage, firstly 
because their reactants are not typically present in large 
quantities at the same times; and secondly, because when 
the reactants H2CO, HCOOH, and the methanol required 
as a precursor for reaction (4) are all abundant, temper- 
atures are much higher (around 100 K), and the associa- 
tion rates are consequently much lower. The second, final, 
peak in the gas-phase methyl formate abundance arises 
from the evaporation of the surface species. 

Figure [3t> plots the grain-surface abundance of 
HCOOCH3. Although not perceptible from the plot, there 
is a subtle shift in the formation route as the methyl 
formate abundance grows. The initial rise at around 
5 — 6 x 10 4 yr is brought about by the surface reaction 
(1) between the radicals HCO and CH3O. This rise takes 
place when CO has strongly evaporated, whilst the eleva- 
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Fig. 3. Chemistry during standard warm-up phase - tf = 2 x 10 5 yr, T — T2(t). Surface species are marked (s). All 
abundances are shown with reference to molecular hydrogen. 



tion of the grain-surface temperature increases the mobil- 
ity of OH and H 2 CO, causing the rate of the reaction 

OH + H 2 CO -» HCO + H 2 

to greatly increase. (The OH radical is produced by the 
photodissociation of water ice, initiated by cosmic ray- 
induced photons). The mobility of the resultant HCO 
has also greatly increased. At these warmer temperatures, 
the hydrogenation of CO has ceased to be important to 
HCO formation, due to the lower abundance of atomic H 
(caused by its faster thermal evaporation). However, when 
CO strongly evaporates at ~ 20 K, the reaction CO + OH 
— > C0 2 + H - the major destruction route of OH - is re- 
moved. This greatly increases the OH abundance, allowing 
much HCO to be formed. At this point in time, although 
atomic H is more scarce than at lower temperatures, the 
reaction 

H + H 2 CO -> CH3O, 



is still the dominant mechanism for CH3O formation. 
H 2 CO is abundant on the grains, having built up during 
the collapse phase. When significant amounts of H 2 CO 
begin to evaporate, at -8 x 10 4 yr (T ~ 40 K), HCO 
and CH3O production is reduced and surface formation of 
HCOOCH3 is inhibited. At this point, reaction (5) starts 
to have a strong effect due to the H 2 CO now present in 
the gas phase, and much of the methyl formate produced 
is then accreted. From this time until HCOOCH3 evap- 
orates, the two processes operate in tandem, with the 
gas-phase/accretion mechanism dominant but gradually 
decreasing in influence. The gas-phase/accretion route is 
responsible for the formation of around 25 % of the total 
HCOOCH3 present on the grain surfaces before evapora- 
tion. 



At early times, HCOOH is formed in the gas phase 
from the dissociative recombination of HCOOH^ , which 
forms by radiative association of HCO + and H 2 0. 
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However, when CH4 evaporates at around 5 — 6 x 10 yr 
(see figure for the gas-phase profile), it can react with 
ions, forming HCOOH^ , which then dissociatively re- 
combines with electrons. This process results in the rise in 
the gas-phase formic acid abundance beginning at 6 x 10 4 
yr. At the same time, as with methyl formate, HCOOH 
surface formation benefits from the increased formation, 
and diffusion rate, of HCO, causing the first rise in surface 
abundance. 

The next rise in the gas-phase HCOOH abundance 
comes about when formaldehyde evaporates, which facil- 
itates more OH formation in the gas phase via reaction 
with atomic oxygen. The OH radical and H2CO then re- 
act in the gas phase to form HCOOH via the process 
OH + H 2 CO -> HCOOH + H. 

Both of these gas-phase mechanisms increase the grain- 
surface HCOOH abundance by subsequent accretion of 
the product. 

On the grain surfaces, reaction (3) becomes most im- 
portant when formaldehyde has evaporated significantly, 
because previously most OH was destroyed by reaction 
with H2CO. After its evaporation, formaldehyde may still 
accrete and then (if it does not immediately evaporate 
again) it may react with OH to form HCO. Sufficient OH 
is still available to further react with HCO via reaction (3) 
because the surface abundance of H2CO is comparatively 
small. At this stage, the dominant OH formation mecha- 
nism is photodissociation of H2O by cosmic ray photons. 
The post-evaporation surface abundance of formaldehyde 
is small compared to its pre-evaporation level, but the 
gas-phase abundance remains high, and so accretion is 
fast and remains steady. Hence, we find that the forma- 
tion of formic acid on surfaces is optimised at a tempera- 
ture where the formaldehyde does not stick to the grains 
(^40 K), but resides in the gas phase and accretes and re- 
evaporates quickly. This allows it to react with OH to form 



enough HCO for reaction (3), whilst not dominating the 
OH chemistry to such an extent that reaction (3) suffers. 
Reaction (3) requires HCO to be formed at a significant 
rate, since HCO quickly evaporates at such temperatures. 
This optimum state for reaction (3) ends when the tem- 
perature increases to the point where OH evaporates (~70 
K). 

However, the grain-surface formation mechanism is 
never dominant; the gas-phase synthesis (OH + H 2 CO) 
is still around ten times faster at this point, due to 
high formaldehyde abundances. The final peak in sur- 
face HCOOH abundance and the small blip before the 
final gas-phase peak are due to increased gas-phase for- 
mation as OH evaporates from the grains. Finally, the 
surface-bound HCOOH evaporates, leaving a high gas- 
phase abundance of ~ 10~ 7 with respect to H2. The ma- 
jority of HCOOH is formed by the gas-phase reaction. 

Dimethyl ether (CH3OCH3) is formed initially on 
grain surfaces by reaction (2). The strong rise up to peak 
surface abundance begins when the surface HNO abun- 
dance falls. This species is until then the main reaction 
partner of CH3, forming CH4 and NO. CH3 is formed on 
the surfaces by the hydrogenation of CH2 following the 
photodissociation of CH4 by cosmic ray-induced photons. 
The fall in HNO allows CH3 abundances to rise, fuelling 
reaction (2). Dimethyl ether continues to be formed by 
the grain-surface mechanism until it evaporates from the 
grains. 

Reaction (2) does not contribute much dimethyl ether 
abundance compared to the final level, and this molecule 
evaporates earlier than methyl formate or formic acid, at 
a time of ~ 10 5 yr. The main formation mechanism is 
therefore the gas-phase route, which consists of the reac- 
tion between methanol and protonated methanol, which is 
fast when methanol finally evaporates, followed by disso- 
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ciative recombination. We use a recombination efficiency 
for producing dimethyl ether of ~1.5 %. 

3.2.2. Other Species 

The gas-phase profile for H2CO (figure^) also shows a 
complicated structure. Again, this is due to the time de- 
pendence of the evaporation of different species. The first 
dip (at ~ 4.5 x 10 4 yr) is due to N2 evaporation, which 
contributes strongly to the NH3 gas-phase abundance; the 
ammonia competes with atomic oxygen to react with CH3 , 
which is required to form H2CO in the gas phase at that 
time. The subsequent evaporations of CO and CH4 then 
push up formaldehyde abundance strongly. The final rise 
in gas-phase abundance takes place as formaldehyde itself 
is desorbed from the grain surfaces. 

The profiles of sulphur-bearing species resulting 
from this treatment show some features of interest. 
Importantly, H2CS is the most abundant sulphur-bearing 
species at late times, in spite of the dominance of H2S on 
the grain surfaces before this. When H 2 S evaporates, it 
is broken down to atomic sulphur by cosmic ray-induced 
photons. Meanwhile, CH3 is abundant, both as a result 
of the high CH4 level in the gas phase, and the cosmic 
ray-induced photodissociation of grain-surface methanol 
followed by the fast evaporation of the resultant CH3. 
Atomic S and CH3 then react in the gas phase to pro- 
duce H2CS. When methanol evaporates, CH3 may still be 
formed by cosmic ray-induced photodissociation, but also 
by dissociative recombination of CHaOHCHj". The final 
H2CS abundance is rather larger than observations might 
suggest, and larger than other models predict. Van der Tak 
et al. (2003) suggest values of 10~ 9 , using radiative trans- 
fer models with observations of high mass star-forming 
regions. Typical H^CSiH^S values are on the order of a 
few. However, that study suggests SO and SO2 should be 



of the same order as H2S and H2CS, which is not the case 
here, except at times of around 8 x 10 4 yr. The observa- 
tional study of IWakelam et all l)2004ri) provides SO, S0 2 
and H2S values on the order of 10~ 6 . 

However, the sulphur chemistry is complex, and is not 
the main focus of this paper. We direct the reader to Viti, 
Garrod & Herbst (2006, in prep.), wh i ch bu ilds on the sul- 
phur chemistry analysis of lViti et alJ (J2QQ4J) , and includes 
grain-surface chemistry. 

3.3. Timescale and Temperature Profile Variations 

In addition to the standard switch-on time of 2 x 10 5 yr, 
we ran additional models with heat-up times of 5 x 10 4 
and 1 x 10 6 yr. Figures 0] and show gas-phase and grain- 
surface abundances of HCOOCH3 and other molecules, 
for the shorter and longer times, respectively. With the 
shorter time to reach a temperature of 200 K, the gas- 
phase routes to HCOOCH3 become dominant over the 
surface reactions, whilst for the longer timescale, the sur- 
face route dominates even more strongly, and also pro- 
duces the most methyl formate of any of the models. 

The major gas-phase route that forms methyl for- 
mate, reaction (5), is strongly dependent on the H 2 CO 
abundance. With the shorter timescale, tf = 5 x 10 4 yr, 
formaldehyde is abundant until the end of the model. This 
is not the case with the longer timescale, tf = 1 x 10 6 yr, 
and here the gas-phase route acts for only a short time 
in comparison to the length of time that the HCOOCH3 
resides on the grain surface, gradually being destroyed by 
cosmic ray-induced photodissociation. However, the pro- 
ductivity of the grain-surface mechanism is enhanced by 
the longer timescale, because there is a longer period be- 
tween the time when HCO becomes both abundant and 
highly mobile (i.e. following CO evaporation), and the 
time when formaldehyde evaporates - the requirements 
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Table 4. Selected gas-phase fractional abundances with respect to H 2 , determined at the time of maximum HCOOCH3 
abundance. 



Species Model 



t/ = 5xl0 4 yr tj = 2 x 10 5 yr t f = 1 x 10 6 yr tj = 2 x 10 5 yr 
T = T 2 (t) T = T 2 (t) T = T 2 {t) T = Tx(t) 



HCOOCH3 (all reactions) 


1.9(- 


-9) 


7.8(- 


-9) 


1.9(- 


-8) 


5.3(-9) 


HCOOCH3 (no gas-phase routes) 


1.0(- 


•10) 


5.5(- 


-9) 


1.9(- 


-8) 


3.2(-9) 


HCOOH 


5.5(- 


-8) 


1.1(- 


-7) 


1.2(- 


-7) 


1.2(-7) 


CH3OCH3 


6.1(- 


-9) 


1.4(- 


-8) 


1.1(- 


-9) 


2.2(-8) 


H 2 CO 


2.3(- 


-5) 


2.3(- 


-7) 


3.9(- 


•10) 


1.5(-8) 


CH3OH 


3.2(- 


-5) 


1.7(- 


-5) 


5.9(- 


-7) 


1.8(-5) 



a a{b) = a x 10 6 



for reaction (1) to be effective. Since the surface formalde- 
hyde abundance is steady until it ultimately evaporates, 
and the OH required to form HCO (via OH + H 2 CO -> 
HCO + H 2 0) is provided steadily by photodissociation of 
water ice, the amount of methyl formate formed by surface 
reactions is directly proportional to the timcscalc. 

For the shorter timescale, the comparatively long life- 
time of H 2 CO in the gas phase means that most HCOOH 
is formed in the gas phase. Hence, shorter timescales 
favour the gas-phase routes, whereas longer timescales 
favour the grain-surface mechanism. 

The final CH3OCH3 abundance suffers with the longer 
timescale, as photodissociation reduces the surface-based 
methanol abundance before it can evaporate and form 
dimethyl ether. However, the early peak in CH3OCH3 
abundance caused by the evaporation of the surface 
species is almost as strong as the later peak induced by 
the gas-phase chemistry. Surface production is increased 
for similar reasons as stated for HCOOH and HCOOCH3. 

We also ran models using the linear temperature de- 
pendence, ATi, shown in figure El with the standard 
switch-on time. The difference between the two profiles 
is that for AT 2 , the increase in temperature is slower at 
early times, but then much faster at later times. Hence, 
with ATi, the hot core spends less time at lower tem- 
peratures and more time at higher temperatures. Because 
the surface reactions (1) - (3) are disfavored for shorter 
timescales, and because these reactions are most active 
from ~20 - 60 K, the ATi profile results in less surface 
formation since it increases more rapidly in this range of 
temperatures. 

4. Discussion 

Our results show that significant abundances of saturated 
organic molecules are produced in hot cores during the 
protostellar warm-up phase. The results also demonstrate 
that the gas-phase and grain-surface chemistry during this 
period are strongly coupled. The large abundances which 
build up on grain surfaces have a large impact on gas- 
phase chemistry when they evaporate. But, also, the prod- 



ucts of that chemistry can re-accrete, changing the sur- 
face chemistry. In addition, the sudden absence of certain 
species from the grains, following evaporation, can have 
a strong effect on the surface chemistry. Since hydrogena- 
tion is inefficient at warmer temperatures, the abundance 
of heavy radicals on the grain surfaces is mainly due to 
the photodissociation of ices such as H 2 0, and to a lesser 
extent CH 3 OH, due to cosmic ray-induced photons. The 
gas-phase and grain-surface chemistries act in different, 
complementary ways, and the strongly coupled system 
yields results which only a full gas-grain treatment can 
produce. 

Table 0] lists the gas-phase abundances obtained with 
various model parameters for methyl formate, formic acid, 
dimethyl ether, formaldehyde, and methanol at the time of 
maximum gaseous methyl formate abundance, which oc- 
curs at or near the end of the warm-up phase. These abun- 
dances derive from both grain-surface and gas-phase reac- 
tions. We have also listed the methyl formate abundances 
obtained without the gas-phase route, since reaction (5) 
is problematic. If we look at the first three columns of the 
table, obtained with varying times for the warm-up phase 
but all with our standard AT 2 profile, we see that dif- 
ferent species show different sensitivities to the timescale, 
which is related to the mass of the protostar. Methyl for- 
mate increases one order of magnitude in abundance from 
the shortest to the longest timescale, with a stronger de- 
pendence if the gas-phase synthesis is removed. On the 
other hand, formic acid has an abundance that does not 
change much. Dimethyl ether shows a rather strange de- 
pendence, first increasing then decreasing as the time in- 
terval increases. The species with the strongest time de- 
pendence is formaldehyde, the abundance of which de- 
creases by almost five orders of magnitude from the short- 
est to the longest time interval considered. The abundance 
of methanol also decreases with increasing time interval, 
although not to as great an extent as that of formaldehyde. 

Let us now consider the agreement between our new 
approach to hot core chemistry and observations of or- 
ganic molecules in these sources. Table |S] shows observa- 
tional results for two hot cores and two hot corinos - IRAS 
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Species 


Orion Compact Ridge 6 


Orion Hot Core 6 


IRAS 16293 c 


IRAS 4A d 


HCOOCH3-A 


3(-8) 


1.4(-8) 


1.7±0.7(-7) 


3.4±1.7(-8) 


HCOOCH3-E 






2.3±0.8(-7) 


3.6±1.7(-8) 


HCOOH 


1.4(-9) 


8(-10) 


~6.2(-8) 


4.6±7.9(-9) 


CH 3 OCH 3 


1.9(-8) 


8(-9) 


2.4±3.7(-7) 


<2.8(-8) 


H 2 CO 


4(-8) 


7(-9) 


6.0(-8) e 


2(-8)' 


CH3OH 


4(-7) 


1.4(-7) 


3.0(-7) e 


<7(-9) 



a a(b) = a x 10 b 

6 Values from &atto^^L (1^9^). Values are given for total HCOOCH3 abundance. 
c Values from Cazaux^j^iL^ (2003). 
d Values from Bottinelli et al. (20041 . 
e Values from Schoier et al. (2002). 
1 Values from Maret et al. (2004) . 



16293 and IRAS 4A - associated with low-mass protostars 
rather than the high-m ass protostars typical of hot cores 
ijBottinelli et alJ l2004'). One can see immediately that ob- 
served abundances for IRAS 16293 arc larger than for the 
other cores. However, there is some disagreement over the 
H 2 column densitie s associated with this source (see e.g. 
IPeeters et aLlEooijl , and this may help to explain the dis- 
crepancy. Indeed, the sizes and morphologies of the emit- 
ting regions of hot cores are not well defined, and therefore 
observationally determined abundances may be strongly 
dependent on beam size. 

For the Orion Compact Ridge and Hot Core, the ob- 
served abundances are best fit by a heat-up interval of 
between 2 x 10 5 yr and 1 x 10 6 yr. Note, though, that 
the abundances listed in table 01 refer to times somewhat 
less than the full interval. The worst disagreement oc- 
curs with formic acid, where the calculated abundance 
is two orders of magnitude too high. One explanation for 
this discrepancy lies in the gas-phase neutral-neutral re- 
action (OH + H 2 CO) leading to HCOOH; there is only 
a small amount of evidence for this channel (see, e.g., 
http://kinetics.nist.gov/index.php). Methanol is best fit 
at time intervals near the upper value while the calcu- 
lated abundances of formaldehyde and dimethyl ether arc 
too low at this time, and methyl formate is not very sen- 
sitive to changes of time in the interval range considered. 
The hot corino IRAS 4A is somewhat richer in methyl 
formate and formic acid than the hot cores, but has less 
methanol. Nevertheless, the level of agreement is not very 
different in the range of heat-up intervals considered. For 
the case of IRAS 16293, the situation is rather unique. 
Here the formic acid is in agreement with our large cal- 
culated abundance but the methyl formate abundance is 
more than an order of magnitude higher than our largest 
value and the dimethyl ether abundance is similarly large. 
The formaldehyde and methanol abundances are not ex- 
ceptionally large, and so rule out a very small heat-up in- 
terval. In other words, both the hot cores and hot corinos 
in our sample are fit best by heat-up intervals characteris- 



tic of low-mass and intermediate-mass stars, although we 
cannot distinguish between the two. 

Our results, while not in excellent agreement with ob- 
servation and while only focusing on abundances of se- 
lected species, add strength to the case that hot cores 
develop over long timescales, spending significant periods 
at low to intermediate temperatures before reaching the 
~100 - 300 K typical of observed objects. Whilst our sug- 
gested heat-up timescales of 2 x 10 5 - 1 x 10 6 yr are quite 
long, our model suggests that observed hot cores could 
have comparatively short post-evaporation ages, where we 
define the "evaporation stage" as that at which CH3OH, 
HCOOH, HCOOCH3, H 2 and NH 3 evaporate, approx- 
imately at the same time (at around 100 K). Note that 
this post-evaporation era is, according to previous models, 
when gas-phase chemistry is supposed t o produce the com- 
plex species considered here (see e.g. Millar ct"aT1 ll99lt 
ICharnlev et~aT]ll992l: ICaselli et alJll993[) . 

Our standard and longer timescale models, tf = 2 x 10 5 
and tf = 1 X 10 6 yr, spend approximately 3 x 10 4 and 
3 x 10 5 yr, respectively, between the evaporation of the 
large molecules and the end of the model. This corre- 
sponds to a change in temperature from ^100 K up to 
200 K. The models largely produce acceptable agreement 
with observed abundances throughout this period, for the 
species which we investigate here, although for the longer 
timescale model methyl formate abundance has begun to 
depreciate seriously after around 10 5 yr post-evaporation. 
The standard model displays acceptable HCOOCH3 levels 
right to the end of the run, corresponding to ^3 x 10 4 yr 
post-evaporation. Our results (see figures El and EJ there- 
fore suggest that hot cores which show HCOOCH3 abun- 
dances of around 10~ 8 with respect to H 2 could have 
post-evaporation ages of anywhere between zero and at 
least 3 x 10 4 yr (in the case of cores associated with 
intermediate- mass protostellar sources), and 1 — 2 x 10 5 yr 
(in the case of cores associated with low-mass protostellar 
sources) . 
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Our model suggests that grain-surface formation of 
HCOOCH3 occurs within a temperature range of ~ 20 
- 40 K. The lower bound in this range corresponds to the 
evaporation of CO, which leaves the OH radical free to 
react with formaldehyde to form HCO. The upper bound 
corresponds to the evaporation of H2CO, which is crucial 
to the formation of both of the radicals required by reac- 
tion (1). Hence, the binding energies of these two species 
are especially important in determining the time period 
available for reaction (1) to take place. The binding en- 
ergy of HCO is much less important, since it is too reactive 
to build up on the grains; its abundance is determined 
by its rate of reaction with other species. Although the 
strength of its diffusion barrier is related to the binding 
energy (Ejj = 0.5Eb), we do not expect the final amount 
of HCOOCH3 produced to be exceptionally sensitive to 
this value. 

Fo llowing; the experi mental evidence o f ICollings et, alJ 
l|2004h for other species, IViti et alJ l|2004 suggested that 
the evaporation of H 2 CO is determined by co-desorption 
with water ice. In this model we allow thermal evaporation 
of H 2 CO according to its desorption energy, less than half 
that of water. However, our model indicates that most 
formaldehyde should be formed at very late times in the 
collapse phase, so that H 2 CO is unlikely to be trapped in 
the water ices below. Using our model to explicitly track 
the deposition of monolayers would allow us to investigate 
this more fully. 

5. Conclusions 

Using a gas-grain network of reactions, we have investi- 
gated the formation of methyl formate in hot cores, along 
with the related species dimethyl ether and formic acid, 
and to a lesser extent formaldehyde and methanol, dur- 
ing the protostellar warm-up and evaporative phase. We 
undertook this work because standard models of hot-core 
chemistry, which rely on gas-phase chemistry at a tem- 
perature of 100 K or so to produce complex molecules 
from precursor methanol, have severe failings (see e.g. 
iHorn et alJl2004|) . The chemistry during the warm-up pe- 
riod occurs both on grain surfaces and in the gas phase. 
We find that the formation of HCOOCH3, by either mech- 
anism, requires a gradual warm-up of temperatures, rather 
than an immediate jump from 10 K to 100 K, as is often 
assumed. In the case of the grain-surface route, reaction 
(1), this is so that radicals may become mobile enough to 
react, without quickly desorbing. Indeed, whilst the time- 
dependent temperature profiles we use are quite specific, 
grain-surface formation of methyl formate should be viable 
as long as temperatures remain in the 20 - 40 K range for a 
long enough period of time. In the case of the only viable 
gas-phase route, reaction (5), it is because the reaction 
rate is strongly temperature dependent - being fastest at 
low temperatures - but requiring the evaporation of H 2 CO 
at intermediate temperatures to become important. The 
plausibility of the gas-phase route for methyl formate iden- 
tified in this work is a matter for further study. 



During the warm-up phase, atomic hydrogen evapo- 
rates very quickly and hydrogenation is much less impor- 
tant than during the cold phase. The production of sur- 
face radicals is largely determined by cosmic ray-induced 
photodissociation of ices. However, in some cases, the low 
surface abundance of atomic H is still sufficient to pro- 
duce radicals, when the other reactant is very abundant, 
e.g. formaldehyde. 

Our results for a variety of species show that both gas- 
phase and grain-surface reactions are strongly coupled. 
Such a strong interaction allows processes to take place 
that could not happen on either the grain surfaces or in the 
gas phase alone. The results of our calculations show that 
towards the end of the warm-up period, large abundances 
of a variety of organic molecules are now present in the gas, 
in agreement with observation. The observed rich, com- 
plex chemistries of hot cores may thus be representative 
of objects in which complex molecules have very recently 
evaporated (at around 100 K), contrary to the predictions 
of models that utilise onl y a gas-phase chemistry to r epro- 
duce all abundances (e.g. iRodgers fc Charnlevl 200lj) . 

The physical model of the warm-up phase, which we 
base on the analysis of Viti et al. (2004, and references 
therein), assumes a finite time period for the central pro- 
tostellar heat source to reach its switch-on time, and there- 
fore for the surrounding gas and dust to vary in tempera- 
ture gradually. This means that the temperature increases 
should propagate outwards, and that the chemistries of the 
inner regions should be further advanced than the outer. 
Therefore, the typical hot core chemistry itself should 
propagate outwards. The radial size of regions which ex- 
hibit significant complex molecule abundances may then 
indicate the speed at which the temperature gradient 
propagates. 

We hope soon to extend our model to include a much 
larger network of surface reactions and species. This set 
will include other astronomically detected molecules such 
as the isomers of methyl formate and dimethyl ether, as 
well as some postulated and (so far) undetected species. 
Surface reaction of heavy radicals will also be extended 
to allow the formation of other large molecules already 
present in the gas-grain code. 
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